Anisotropic Fermi superfluid via p-wave Feshbach resonance 



in . 
o ■ 
o ■ 

<n : 

Or 
<■ 



5-1 

43 



O 
O 



> 

oo 

(N 
O 

in 
o 
-i— > 



T3 

C 

O 

o 



X 



Chi-Ho Cheng and S.-K. Yip 

Institute of Physics, Academia Sinica, Nankang, Taipei 115, Taiwan 
(Dated: February 2, 2008) 

We investigate theoretically Fermionic superfluidity induced by Feshbach resonance in the orbital 
p-wave channel. We show that, due to the dipole interaction, the pairing is extremely anisotropic. 
When this dipole interaction is relatively strong, the pairing has symmetry k z . When it is relatively 
weak, it is of symmetry k z +if3k y (up to a rotation about z, here f3 < 1). A phase transition between 
these two states can occur under a change in the magnetic field or the density of the gas. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.90.+q 



In the earlier investigations of trapped Bose and Fermi 
gases, the interactions between particles are in general 
weak since the gases are typically dilute, with the inter- 
particle distances much larger than the scattering a be- 
tween the particles (see, e.g., |l|,|2j and reference therein). 
Recently however, it has been recognized that the in- 
teraction among particles can be tuned with the use of 
Feshbach resonances [3J. In the case of two Fermion 
species and a Feshbach resonance between them, one 
can tune, by varying the magnetic field and hence the 
resonance level energy, the interaction from (its "back- 
ground" value through) weakly to strongly attractive. 
The ground state of the system is expected to evolve from 
a Bardeen-Cooper-Schrieffer (BCS) superfluid with long- 
range (compared with interparticle distances) Cooper 
pairing to a Bose-Einstein condensate (BEC) of tightly 
bound molecules. £| This cross-over has already been 
a subject of intensive and experimental investigations. 



Practically all the above investigations are for the case 
where the resonance involves an s-wave bound state in 
the "closed" channel. Thus the aforementioned states 
involve pairing only in the s-wave state. More recently, it 
has been demonstrated experimentally also that p-wave 
Feshbach resonances exist 0, S • This thus raises the 
possibility of superfluids with p-wave Cooper pairing or 
BEC of p-wave molecules [inj . 

Consider now therefore a single Fermion species with 
an effective interaction induced by the p-wave Feshbach 
resonance. This brings in mind the well-known system of 
superfluid 3 He. [ll( 3 He has (nuclear) spin 1/2, in gen- 
eral not spin-polarized and basically spatially isotropic. 
Back in the 60's, Anderson and Morel ^3 investigated 
theoretically the superfluid state for this system by as- 
suming that the pairing exists only between the same 
(say |) species. They found that the ground state of this 
system corresponds to Cooper pairing in the I = 1, m = 1 
channel, i.e., the pair-wavefunction has the symmetry of 
the spherical harmonics Y 1 (k) cx [k x + ik y ) (or its spatial 
rotation) . This state is realized in the 3 He A-phase and is 
known as the "axial" phase One can thus conclude 
immediately that for a spin-polarized but otherwise spa- 
tially isotropic system, the pairing is again expected to 
be in the Yi(k) state. This argument is in agreement 



with the findings of Ref 10]. 

However, in our atomic system of interest, the interac- 
tion is far from isotropic. In particular, as demonstrated 
and explained by Ticknor et al || for the case of 40 K, 
due to the magnetic dipole of the alkali atoms, the 1 = 1, 
to = resonance occurs at a higher magnetic field than 
the m = ±1 ones. For a given magnetic field, the in- 
duced effective interaction is actually more attractive in 
the 777 = channel than the m = ±1 ones. It is now thus 
non-trivial what the symmetry of ground state should be. 
This is the specific question on which we would concen- 
trate in this paper. 

Our findings are basically as follows. For a sufficiently 
dilute Fermi gas or for the case where the m = and 
777 = ±1 resonances are sufficiently far apart, the pairing 
occurs only in the m = channel. That is, the ground 
state is a BCS state with Cooper pairing of symmetry 
Y®(k) cx k z for "large" magnetic fields and a BEC state 
of 777 = Bosons for magnetic fields sufficiently below the 
resonance(s). There is no pairing or Boson in the m = ±1 
channel. The orbital symmetry of this pairing is the same 
as the "polar" phase in the 3 He literature For a 

sufficiently dense gas or for the case where the resonances 
are sufficiently closed to each other, the BCS case at high 
magnetic field corresponds to a state with Cooper pairing 



i/3(k ■ a) 



Here z is along the magnetic field 
direction and a is a vector in the plane perpendicular to z, 
and (3 is a number less than unity and is field dependent. 
This state is thus intermediate between the polar and 
axial phases. At lower fields, (3 may vanish and again the 
BCS pairing or the BEC condensation is again entirely 
in the m — channel. These results are summarized in 
Figs. HJandEl 

We begin with the Hamiltonian H = Hf + Hf, + H a 
where 
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This Hamiltonian is a generalization of the ones al- 
ready commonly employed 0, Q, Hi! flfij ] for s-wave Fes- 
hbach resonances to the case of several I = 1, m — 0, ±1 
closed channels. Hf and Hb are the Hamiltonians for the 
free Fermions and Bosons (particle pairs in the "closed" 
channels) respectively, and H a represents the Feshbach 
coupling. et£ is the annihilation operator for a Fermion 
with momentum k, the corresponding operator for a 
Boson with angular momentum I = 1 and z-axis projec- 
tion m with momentum q, [i and M are the chemical 
potential and mass of the Fermions. In H a , the factor 
Y 1 m (fc) reflects the symmetry of the I, m bound state and 
the linear factor in k arises from the small momentum ap- 
proximation for the coupling. S m is the (bare) detuning 
of the energy of the closed channel with angular momen- 
tum projection m, and a m the corresponding coupling 
constant. For low energy scattering between a pair of 
particles of momentum ±fc, the scattering amplitude in 
the I = 1, m partial wave can be parameterized as, using 
the same notations as ref Ig, 
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v m has the dimension of a volume and c m an inverse 
length. The magnetic field dependent parameters v m , 
c m are in principle available experimentally and have al- 
ready been measured for K for atoms in the |/, mt >= 
|9/2, — 7/2 > hyperfine states, [8j The bare parameters 
S m and a m can be related to the physical parameters v m , 
c m by considering scattering between two Fermions using 
the Hamiltonian given in eqs (HEJ). 1JJ These relations 
are: 
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momenta. Below we shall express all physical quantities 
in v m and c m in the final expressions and omit these 
explicit cutoffs. 

With Feshbach resonance for the sub-channel m, X/v m 
is field dependent, vanishing at a field B^. In contrast, 
c m has a definite sign. For the ease of discussions, we 
shall assume that c m < and field independent, —l/v m 
is an increasing function of field —l/v m > (<)0 for B > 
(<)B* n , as in the case of 40 K. (This corresponds to the 
case where S m is an increasing function of field and a m 
weakly field dependent, c.f. eqs (JSJ and (JSJl). For B < 
B^, a bound state appears. The energy of this bound 
state is given by —tb.m = —h 2 K m 2 /M with k = lK m 
being a pole for f m (k). For small detuning below the 
resonance, k 2 ^ = l/[(— c m )(v m )]. Since l/v m should be 
roughly linear in B near the resonance, tb,m increases 
linearly with (B^ — B) (in contrast to s-wave, where it 
is quadratic). 

Moreover, as explained in due to the dipole interac- 
tion, Bq > B{ = B*_ x . Thus in the field range of interest, 
—l/v\ — —1/v-i > —1/vo- We can say that, at a given 
field, the effective interaction between the Fermions is 
less attractive for relative angular momentum projections 
m = ±1 than m = 0. 

Now we proceed to find the ground state for the many- 
body problem. We assume a mean-field theory and re- 
place bfrm by c-numbers. Only its q — value is non- 
vanishing. It is convenient to introduce the symbols 
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so that H a becomes Aia_gog + /i.e.. The Fermionic 
part of the Hamiltonian can be solved by Bogoliubov 
transformation. The value for &o,m can also be easily 
found since the Hamiltonian is quadratic in this variable. 
In terms of D m , we get 



Here = j^j. The divergent sums over k on the right- 
hand-sides of the above two equations can be regulated 
either by introducing a cut-off or invoking the fact that 
the coupling a rn must actually decay to zero at large 
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Using eq JSJ) and I© we obtain 
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and a corresponding equation with m = 1 m = — 1. 
Here 
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These equations are to be solved together with the num- 
ber equation 
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(at 0jE ) is given by v\ = |(1 - 

from eq JSJ is linear in fc, the sum in eq (| 1 31) over 

k is formally divergent due to the large k contributions. 
However, we can regularize it by employing again eq ffijl. 
We have finally 
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Eq (|10|) . (flip) . (|14fl are our principal equations, with 
parameters characterizing the Feshbach resonances ex- 
pressed entirely in v m and c m . These equations deter- 
mine the order parameters D m and chemical potential fi 
for given density n and " interaction parameters" v m and 
c m . For simplicity, in writing these equations we have 
already dropped the terms with explicit l/|a m | 2 factors. 
These terms are small under the " wide-resonance" regime 

0(b), HQ. 

Since the interaction is less attractive for angular mo- 
mentum projections m = ±1, for sufficiently large differ- 
ence between —1/vq and — l/u±i we expect (and verify 
below) that the pairing is entirely in the m = par- 
tial wave. We thus first begin our analysis by assuming 
that only Do is non-vanishing, eqs IjlUI) and (|14fl can be 
solved simultaneously similar to the s-wave case. It is 
convenient to express the results in dimensionless form. 
We define jl = /x/e F , D m = D m /v F , c m = n~ 1/3 c m , 
and v m = nv m where cf = ft 2 fcp/2M, vf = k^/M, and 
/cf = 67r 2 n. The results are as shown in Fig. ^ (for 
the case Co = c\ = —100, see below for the reason of 
this choice). In the BCS regime ( large —1/vq ^s> 1 



or B — £>q, not shown explicitly), p, — > 1, and Do, be- 
ing proportional to the magnitude of the BCS gap, is 
« 1. In the BEC ( — 1/io < -1 or large and negative 
B — Bq), /i is approximately — e^/2 and Do approaches a 
constant. This latter value can be obtained from eq (|14|> 
as (32V3) 1/3 (-5 r 1/2 - 

The "cross-over" behavior in Fig ^ is analogous to 
the s-wave case, where the corresponding x-axis is x — 

— l/(n 1//3 a) where a is the s-wave scattering length. Note 
here Do has the dimension of (energy x inverse length) 
and behaves differently from the s-wave A in the BEC 
limit. We have also performed calculations for other val- 
ues of Co . The size of the crossover region is roughly pro- 
portional to the value of Co- For example, for cq = —200, 
the corresponding results can be captured well by replac- 
ing the x-axis by — 1/2% and dividing Do by 1/ y2 in Fig 
^ The result in Figgis similar to that in Ref |l0|, even 
though the latter actually studied a different (D\ ^= 0, 
Do = D-i = 0) state. 

The above behavior applies only to sufficiently large 

— l/v±i — (— 1/uo) > 0. When this difference is suffi- 
ciently small, D±i will become finite. The critical value 
for l/v±i, denoted by l/t>±i, can be found by putting 
D±i — in eq iJTUjl and linearizing eq ifTTJl in D±\. Wc 
obtain, for ci = Co, 
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In the BCS limit (— l/v >• 1), the first term in eq 1(15)1 
is negligible whereas the second term becomes a constant 
independent of Do- From this we get —1/v* + 1/vq — » 
12tt w 37.7. In BEC limit (-l/v < -1), the main 
contribution comes from the first term in eq (|15|l . Us- 
ing the aforementioned asymptotic values of Do we get 
-1/5J + l/v -> 48tt/5 w 30.2. -l/v\ + l/v is shown 
as the thick black line (Do = 0) in Fig [21 

For —1/vi + 1/io less than this critical value, D±\ are 
finite. Eqs l|10[). (jl 1|> involve three complex unknowns 
D m . By gauge invariance we can always choose Do to 
be real. Under this choice, the solutions we found be- 
long to the class D 1 = D*_ v Writing D x = \Di\e ix , A £ 
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then has the angular dependence oc D k z + iy/2\Di\k ■ a 
oc {k z + ij3k ■ a) where a = (cosx)y + (sinx)i is a unit 
vector perpendicular to z and (3 = \/2\Di/D \. A partic- 
ular solution is given by the case where Di and D_i are 
both real where a = y. The other solutions are simply 
related to this one by a rotation about z. Without loss 
of generality we shall therefore pretend that D m 's are all 
real. 

The contours of the order parameters D±\ are also 
shown in Fig [5] In the D±i = phase, the state is rota- 
tionally invariant about z, whereas this symmetry is bro- 
ken in the D±i ^ phase. There is a (quantum) phase 
transition between these two phases when one crosses the 
line -1/v* + 1/vq. 

For the 40 K case studied in |8j , the Feshbach resonances 
are at B% w 198.8G and B{ w 198.4G. There, a and 
Co are both only weakly field dependent and are approxi- 
mately given by — 0.02a " 1 . Our choice of c = —100 above 
corresponds to a density of roughly 6.7 x 10 13 cm~ 3 . Near 
the resonant fields, —\/v\ + l/v ~ 2.1 x 10~ 8 a,7 3 and is 
roughly field independent. Thus the density determines 
the values for both cq,i and — l/«i + 1/uo while varying 
the magnetic field corresponds roughly to moving along 
a horizontal line on our phase diagram of Fig. [3 (with 
increasing field towards the right and the distance of the 
line from the x-axis proportional to n^ 1 ). 

While preparing this manuscript, we become aware 
of Ref (jjj which studies essentially the same problem. 
Whereas our results agree with |17| for very large and 
very small splitting between the resonances, the conclu- 
sions differ in the intermediate splitting regime. Our pre- 
diction is that the state is ~ k z + i(3k y on the BCS side 
whereas it should be k z on the BEC side. Their con- 
clusion is the opposite. The reason for this disagreement 
is not yet understood. We believe that our results are 
more reasonable. For large positive detuning, the split- 
ting should be less relevant and the pairing state should 
resemble more that of the isotropic system. On the BEC 
side, the system should be closer to a Bose condensate of 
lowest energy (k z ) molecules. 

In conclusion, we have shown that p-wave Feshbach 
resonance in general leads to anisotropic Fermi superflu- 
ids. The symmetry of the ground state depends on both 
the density and magnetic fields. 
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FIG. 1: The dimensionless parameters Do and fi as functions 
of — 1/vq. Co = Si = —100 and D±i = in this case. The 
dashed line represents — £(,/2eF. 
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FIG. 2: Contour plot of D±i as a function of 
— 1/vo for Co = 5i = —100. The line for Z?±i 
to the critical value —1/v* + 1/vq. 
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